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Q^l We present the broad band UBVI CCD photometric investigations in 

P ■ the region of the two open clusters Haffner 11 and Czernik 31. The radii of 

j^ . the clusters are determined as 3'. 5 and 3'.0 for Haffner 11 and Czernik 31 

\ respectively. Using two colour {U — B) versus {B ~V) diagram we determine 

^ ■ the reddening E(B — V) = 0.50±0.05 mag and 0.48±0.05 mag for the cluster 

^ ■ Haffner 11 and Czernik 31 respectively. Using 2MASS JHKg and optical 

"^ : data, we determined E{J-K) = 0.27±0.06 mag and E{V-K) = 1.37±0.06 

O ■ for Haffner 11 and E{J - K) = 0.26±0.08 mag and Eiy - K) = 1.32±0.08 

mag for Czernik 31. Our analysis indicate normal interstellar extinction law 
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i^ , in the direction of both the clusters. Distance of the clusters is determined 

j^ \ as 5.8±0.5 Kpc for Haffner 11 and 3.2±0.3 Kpc for Czernik 31 by comparing 

the ZAMS with the CM diagram of the clusters. The age of the cluster has 

been estimated as 800±100 Myr for Haffner 11 and 160±40 Myr for Czernik 

31 using the stellar isochrones of metallicity Z = 0.019. 
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1. Introduction 

Open clusters (OCLs) are located in the Galactic disc. It is important to 
determine their properties and spatial distribution to understand the struc- 
ture and evolution of Milky Way. Therefore, investigation of the OCLs and, 
in particular, the estimation of their physical parameters like age, distance, 
reddening, size and metallicity are very valuable. 

Considering the above we decided to study two open star clusters namely 
Haffner 11 (Ha 11) and Czernik 31 (Cz 31). Ha 11 was first discovered by 
Haffner (1957). Further, van den Bergh et al. (1975) studied this cluster pho- 
tographically and derived a diameter of 5' and richness of type M(moderate). 
Cz 31 was discovered by Czernik (1966) using the charts of Palomar sky atlas. 
This cluster is classified as a Trumpler class IV2p by Ruprecht (1966). 

The main purpose of the present investigation is to estimate the funda- 
mental parameters of the clusters Ha 11 and Cz 31 using UBVI CCD data. 
In the literature, both the clusters are poorly studied. The basic parameters 
taken from WEBDA and Bias et al. (2002) are listed in Table [H Both the 
clusters are located in the third Galactic quadrant and towards the Galactic 
anticenter direction. These objects are very useful to understand the disc 
subsystem to which the cluster belongs. Table [1] shows that the clusters 
are of intermediate and young age objects which can be used to probe the 
chemical evolution and star formation history of the disc. 



In the present study, we provide UBVI CCD photometry for Ha 11 and 
Cz 31 and determine their basic parameters using optical and near-IR 2MASS 
JHKs data. CCD optical data used in the present analysis is obtained for 
the first time for Cz 31. 

2. Observations and data reduction 

The CCD broad band UBVI images were collected using 2Kx2K CCD 
system at the //13 Cassegrain focus of the Sumpurnanand 104-cm telescope 
located at ARIES, Manora peak, Nainital, India. The used CCD has 24 /xm 
square pixel size, resulting in a scale of 0".36 pixel"^ and a square field of 
view of 12'. 6 size. The CCD gain was 10 e^/ADU while the read out noise 
was 5.3 e^. In order to improve the S/N ratio, the observations were taken 
in 2x2 pixel binning mode. The observations include several exposures in 
each filters. The details of observations are listed in Table [2l Dates of 
observations together with the filters used and the corresponding exposure 
time are provided in Table [21 Identification maps of the observed regions of 
the clusters are shown in Fig[Tl 

A number of bias and flat-field frames were taken during the observa- 
tions. Flat-field exposures were made of the twilight sky in each filter. Cor- 
rections for bias and flat-field were performed using the standard IRAFcl 
procedure. The subsequent data reduction and analysis were done using 
the DAOPHOT software (Stetson (1987), (1992)). The stellar photometric 
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routine of DAOPHOT was used for the magnitude determination. These 
magnitudes were cahbrated using Landolt (1992) standards. The details of 
data reductions are described in the previous papers (Pandey et al. (1997); 
Durgapal et al. (1997, 2001)). The instrumental magnitudes were obtained 
through quadratically varying Point Spread Function (PSF). For each filter, 
the stars have been aligned to that of a reference frame, the deepest one. 
The average instrumental magnitude was derived in each filter. 

We observed the standard field PG 0942+051 and SA 98 (Landolt 1992) 
in U, B, V and I for calibrating the observations of Ha 11 and Cz 31. The 
standard stars used in the calibrations have brightness and colours range 
11.93 < l^ < 15.67 and -0.29 < (B - V) < 2.19 respectively, thus covering 
the bulk of the cluster stars. For the extinction coefficients we assumed the 
typical values for the ARIES site. 

For translating the instrumental magnitude to the standard magnitude, 
the calibration equations are as follows: 

u = U + Zu + Cu{U -B) + kuX 

b = B + ZB + Cb{B -V) + ksX 

v = V + Zv + Cv{B -V) + kyX 

i = I + Zj + Ci{V -I) + kiX 

where u, b, v and i are the aperture instrumental magnitudes and U, B, V 
and I are the standard magnitudes and X is the airmass. The color coeffi- 
cients (C) and zeropoints (Z) for different filters are listed in Table O The 
errors in zero points and colour coefficients are ~ 0.01 mag for B, V and / fil- 
ters. The internal errors in magnitude derived from DAOPHOT are plotted 



against V magnitude in Fig. HJ This figure shows that photometric error is 
< 0.01 mag at y ~ 19*^ mag for B, V and I filters while error is < 0.05 mag 
for U filter at V ~ 17*'' mag. Photometric global (DAOPHOT+Calibration) 
errors have also been determined by following Patat et al. (2001) and are 
listed in Table Si For the V filter, the errors are 0.05 and 0.06 at V ^ 
17.0 and 20.0 mag, respectively. The final photometric data are available in 
electronic form at the WEBDA site g and also from the authors. 

The X and Y coordinates of the stars in the observed region of the clus- 
ters have been converted to right ascension (RA) and declination (DEC) of 
J2000. In order to obtain an astrometric solution we use the SkyCat tool and 
Guide Star Catalogue v2 (GSC-2) at the European Southern Observatory. 
This way we considered 100 bright stars per field for which we have both 
celestial coordinates on the GSC-2 and the corresponding pixel coordinates. 
By using CCMAP and CCTRAN in IRAF, we estimated the transformation 
and compute the individual celestial coordinates for all the detected stars. 
The transformation have an rms value of about 0.10 arcsec for RA and DEC. 

3. Analysis of the data 

3.1. Cluster radius and radial stellar surface density 

To determine the radius of the cluster, we derive the surface stellar density 
by performing star counts in concentric rings around the cluster centre listed 
in Table m and then divided by their respective areas. In Fig. Owe show the 
density profile for the cluster Ha 11 and Cz 31. The radial density profile 
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of the cluster Ha 11 shows a flattening around r ~ 3.5 arcmin and started 
merging with background stellar density. Therefore, we consider 3.5 arcmin 
as the cluster radius. Our estimate of radius is more than the value listed 
in Table [1] The radial density profile for Cz 31 is decreasing smoothly and 
~ 3.0 arcmin, it is mixing with field stars density. We consider 3.0 arcmin 
radius for Cz 31 which is larger than the value listed in Table [1] 

3.2. Colour-magnitude diagrams 

The V, {B — V) CMDs of the cluster and field region is shown in Fig. H] 
for the cluster Ha 11 and Cz 31. Stars falling within the cluster radius are 
considered as cluster region stars while those outside the radius are assumed 
as field region stars. To get the clear features in the CMD, we consider the 
stars as cluster member within cluster radius. Field region CMDs shown in 
Fig. m are clearly dominated by foreground/background stars. 

Haffner 11: The V, {B — V) CMD of this cluster show a main-sequence 
(MS) extending from V ^ IQ mag, where the turn-off is located, down to 
\^ ~ 18 mag. After V ~ 18 mag, field stars are clearly dominating and MS is 
merging in the field stars. Few red giant stars are also visible around l^ ~ 16 
and {B — V) ^ 1.5 mag in the CMD. The morphology of the CMD indicate 
that it is a typical intermediate age open cluster. 

Czernik 31: A close inspection of the CMD exhibits a poorly populated 
MS extending from V^ ~ 12 mag down to V ~ 17 mag. The MS fainter 
than V ^ 11 mag has more scatter and field star contamination is also more 
evident. Because of that it is hard to separate the cluster members from 
field stars only on the basis of closeness to the main populated area of the 
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CMD. The CMD of this cluster looks like a poorly populated young open 
star cluster. 

3.3. Colour-colour diagram 

To determine the reddening of the clusters we plot {U — B) versus {B — V) 
colour-colour diagrams in Fig |5] using probable members of the cluster. The 
intrinsic zero-age main-sequence (ZAMS) given by Schmidt-Kaler (1982) is 
shown by the solid line, whereas the dashed line is the same ZAMS which 
is shifted by assuming the slope of reddening E{U — B)/E{B — V) as 0.72. 
The ZAMS shifted to the MS stars provides a mean value oi E{B — V) = 
0.50 ±0.05 mag for Ha 11 and E{B-V) = 0.48 ±0.05 for Cz 31. Our derived 
values of reddening is in agreement with the value listed in Table [T] for Ha 
11. Bica et al. (2005) estimated the reddening E{B - V) = 0.06 ± 0.03 for 
Cz 31 using 2MASS JHKg photometry. Our derived value of reddening is 
much higher than Bica et al. (2005). 

3.4. Interstellar extinction in near-IR 

The near-IR data is available for both the clusters in 2MASS catalogue. 
The JHKg data in combination with optical data has been used to study the 
interstellar extinction. The Kg magnitudes are converted into K magnitude 
following Persson et al. (1998). The (J — K) versus {V — K) diagram for 
the clusters under study is shown in Fig. O The ZAMS shown by the solid 
line is taken from Caldwell et al. (1993). The same ZAMS is shifted by 
E{J - K) = 0.27 ± 0.06 mag and E{V - K) = 1.37 ± 0.06 mag for Ha 11 
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and E{J-K) = 0.26 ±0.08 mag and ^(^-7^) = 1.32 ±0.08 mag for Cz 31 
and shown by dotted line. The ratios ^L~^l ~ 0.20 ± 0.10 for both Ha 11 
and Cz 31 is in good agreement with the normal interstellar extinction value 
of 0.19 suggested by Cardelli et al. (1989). However, scattering is larger due 
to error in JHK data. 

Using the relation R = 1.1E{V - K)/E{B - V) given by Whittet et 
al. (1980) we studied about the nature of interstellar extinction law in the 
direction of both the clusters. The values of R derived in this way are 3.0 for 
both the cluster Ha 11 and Cz 31. Based on this analysis, we can conclude 
that interstellar extinction law is normal in the direction of both the clusters. 

3.5. Age and distance to the clusters 

The age of a cluster is determine by comparing the theoretical stellar 
evolutionary isochrones given by Girardi et al. (2000) for Z = 0.019 with 
the observed CMDs of the clusters shown in Fig. [71 The detailed shape 
and position of the different features in the CMD depend mostly on age and 
metallicity and also on reddening and distance. 

We surveyed different age and metallicity isochrones to get the best fit 
isochrones to all the CMD features in V, (U-B); V, {B - V) and V, {V - 1) 
diagrams. To get the clear features in the CMDs, we consider only those stars 
which lie within the cluster radius as derived in Section I3TT1 The uncertainties 
in age and distance reflect the range that allows a reasonable fit to the cluster 
CMD. 

Haffner 11: In Fig. [3, we show the fitting of isochrones to V, {U — 
B); V, {B — V) and V, (V — /) CMDs using the reddening derived in Section 
We plotted only those stars which lie within the cluster radius derived 
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in Section IXn The isochrones of different age (log(age)=8.85, 8.90 and 8.95) 
and Z = 0.019 have been superimposed on the CMDs. The overall fit is good 
for log(age)=8.90 (middle isochrone). The detailed shape of the MS, TO and 
giant phase are reproduced. The best fitting isochrone provides an age of 
800 ± 100 Myr. The inferred distance modulus (m — M) = 15.36 ± 0.2 mag 
provide a heliocentric distance 5.8 ± 0.5 kpc. The Galactocentric distance is 
12.3 kpc, which is determined by assuming 8.5 kpc as the distance of the Sun 
to the Galactic center. The distance determined in the present study is in 
agreement with the value listed in Table 1. The Galactocentric coordinates 
are X = -5.1 kpc, Y = 2.6 kpc and Z = -0.36 kpc. Using 2MASS JHR^ 
photometry, Bica et al. (2005) derived a heliocentric distance 5.2±0.2 kpc 
and age 890±150 Myr, which are in agreement with the present estimate of 
distance and age. 

Czernik 31: In Fig. [71 we superimpose isochrones with Z = 0.019 and 
for 3 ages (log(age)=8.10, 8.20 and 8.30) in V, {U - B);V, {B - V) and 
V, {y — I) CMDs. The best fit isochrone is found for an age of 160 ± 50 Myr, 
where the associated error has been derived by trying different age isochrones. 
On average, we obtained distance modulus {m — M) = 14.0 ± 0.2 mag. The 
errors have been derived by displacing the best fit isochrone back and forth in 
the distance modulus direction and exploring the acceptable value of distance 
modulus. The estimated distance modulus provides a heliocentric distance 
3.2 ± 0.3 kpc. The Galactocentric coordinates are X = —2.7 kpc, Y = 1.8 
kpc and Z = 0.01 kpc. The Galactocentric distance of the cluster is 10.6 kpc 
towards the anti-center direction. 



3.6. Near-IR colour-magnitude diagrams 

Using optical and near-IR data we redetermined distance and age of the 
clusters under study. We plot V, {V - K) and K, (J - K) CMDs for Ha 11 
and Cz 31 in Fig. [HI The theoretical isochrone given by Girardi et al. (2000) 
for Z = 0.019 have been overplotted for log(age)=8.9 and 8.1 in the CMDs of 
Ha 11 and Cz 31 respectively. The apparent distance moduli {m — M)v,{v-K) 
and {m — M)j j^j^k) turn out to be 15.4 ± 0.3 and 13.7 ± 0.3 for Ha 11 and 
14.0 ± 0.3 and 12.6 ± 0.3 for Cz 31. Using the reddening derived in Section 
13. 4[ we derive a distance of 5.4 ± 0.5 and 3.3 ± 0.3 kpc for Ha 11 and Cz 31 
respectively. Both, age and distance determination for the present clusters 
is in agreement with the estimates presented in Section 13.51 

4. Conclusions 

We studied two open star clusters Haffner 11 and Czernik 31 using UBVI 
CCD optical and 2MASS JHKg data. The results are summarized in Table 
13 The main findings of our analysis are given below. 

1. The radii of the clusters are obtained as 3'. 5 and 3'.0 which corresponds 
to 5.9 pc and 2.8 pc, respectively, at a distance of the cluster for Ha 11 
and Cz 31. 

2. From the two colour {U — B) versus {B — V) diagram, we estimated 
E{B -V) = 0.50 ± 0.05 mag for Ha 11 and 0.48 ± 0.05 mag for Cz 31. 

3. Distance to the cluster Ha 11 and Cz 31 are determined as 5.8±0.5 and 
3.2±0.3 Kpc respectively. Age of 800±100 Myr and 160±40 Myr are 
determined for the clusters Ha 11 and Cz 31 respectively by comparing 
the isochrones oi Z = 0.019 given by Girardi et al. (2000). These 
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distances and ages are supported by the values derived by combining 
optical and near-IR data. 
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Figure 1: Finding chart of the stars in the cluster and field regions of Haffner 11 and 
Czernik 31. The circle represents the cluster size. 



13 




10 12 14 16 18 20 



Figure 2: Photometric errors in U, B,V and / magnitudes against V magnitude. 
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Figure 3: Surface density distribution of stars in the field of the eluster Haffncr 11 and 
Czernik 31. Errors are determined from sampling statistics(=^= where N is the number 
of stars used in the density estimation at that point). Dotted line represents the level of 
field stars density. 
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Figure 4: The V, {B — V) CMD for the cluster Haffner 11 and Czernik 31 using stars 
within the cluster radius. Stars outside the cluster radius are also plotted as field region 
stars. 
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Figure 5: The ([/ — B) versus (B — V) colour-colour diagram for the clusters under 
study. The continuous curve represents locus of Schmidt-Kaler's (1982) ZAMS for solar 
metallicity. The dashed lines are the same ZAMS shifted by the values given in the text. 
The solid arrow indicates the reddening vector. 
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Figure 6: The plot of (J — K) versus (F — K) colour-colour diagram of the clusters for 
the stars within the cluster radius. The solid line is the ZAMS taken from Caldwell et al. 
(1993). The dotted lines are the ZAMS shifted by the values given in the text. 
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Figure 7: The V, (U-B); V, {B - V) and V, {V - I) CMDs for the cluster Haffner 
11 and Czernik 31 usmg stars within cluster radius. The different lines are the differ- 
ent age isochrones taken from Girardi et al. (2000). Three isochroncs of different age 
(log(age)=8.85, 8.90 and 8.95) of metalhcity Z = 0.019 are shown in the CMD of Haffner 
11. The isochrones of log(age)=8.10, 8.20 and 8.30 and Z = 0.019 are plotted in the CMD 
of Czernik 31. 
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Figure 8: The V versus {V — K) and K versus (J — K) CM diagram of the clusters using 
probable cluster members. The solid curve represent the isochrones of log(age) = 8.90 for 
Haffner 11 and log(age) — 8.20 for Czernik 31 taken from Girardi et al. (2000) for Z = 
0.019. 
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Table 1: Fundamental parameters of the clusters taken from WEBDA and Dias et al. 
(2002). 



Name 


tt2000 


<52000 


/ 


b 


Dia 


E{B-V) 


D 


log(age) 




h:m:s 


d:m:s 


(deg) 


(deg) 


w 


(mag) 


(pc) 




Haffner 11 


07:36:25 


-27:43:00 


242.4 


-3.5 


5 


0.57 


6100 


8.70 


Czernik 31 


07:36:59 


-20:32:00 


236.3 


0.3 


5 


0.06 


2200 


8.25 
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Table 2: Log of observations, with dates and exposure times for each passband. 



Band Exposure Time Date 

(in seconds) 
Haffner 11 
U 1500x2, 300x2 6*^ Feb 2010 
B 1200x2,240x1 

V 900x2, 180x2 
/ 300x2,60x1 

Czernik 31 
U 1200x1, 240x2 5^'' March 2010 
B 900x1, 180x2 

V 600x1, 120x2 
/ 300x2, 120x1 



22 



Table 3: Derived Standardization coefficients and its errors. 



Filter 


Colour Coeff. (C) 


Zeropoint {Z) 




Haffner 11 




U 


-0.040±0.09 


7.35±0.05 


B 


+0.001±0.01 


5.29±0.01 


V 


-0.104±0.01 


4.99±0.01 


I 


-0.125±0.01 
Czernik 31 


5.43±0.01 


U 


-0.03±0.02 


7.85±0.01 


B 


-0.03±0.01 


5.64±0.01 


V 


-0.06±0.01 


5.21±0.01 


I 


-0.06±0.01 


5.48±0.01 
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Table 4: The rms global photometric errors as a function of magnitude. 



V ay (Tb cfi (^u 

9-11 0.03 0.05 0.03 0.05 

11-12 0.05 0.05 0.03 0.05 

12-13 0.05 0.05 0.05 0.06 

13-14 0.05 0.05 0.05 0.06 

14-15 0.05 0.05 0.05 0.06 

15-16 0.05 0.05 0.05 0.08 

16-17 0.05 0.05 0.05 0.11 

17-18 0.06 0.06 0.06 0.13 
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Table 5: Derived fundamental parameters of the clusters under study. Rgc is the Galac- 
tocentric distance while X, Y and Z are the Galactocentric coordinates of the clusters. 
The coordinate system is such that the Y-axis connects the Sun to the Galactic Centre, 
while the X-axis is perpendicular to that. Y is positive towards the Galactic anticentre, 
and X is positive in the first and second Galactic quadrants (Lynga 1982). 



Name 


Radius 


E{B - V) 


Distance 


X 


Y 


Z 


Rgc 


Age 




(arcmin) 


(mag) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(Myr) 


Haffner 11 


3.5 


0.50 ± 0.05 


5.8 ±0.5 


-5.1 


2.6 


-0.36 


12.3 


800 ± 100 


Czernik 31 


3.0 


0.48 ± 0.05 


3.2 ±0.3 


-2.7 


1.8 


0.01 


10.6 


160 ± 40 
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